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Abstract 



^ ■ We study, using dissipative particle dynamics simulations, the effect of active lipid flip-flop on 

(N 

0^ . model fluid bilayer membranes. We consider both cases of symmetric as well as asymmetric flip- 

o 

•^ \ flops. Symmetric flip-flop leads to a steady state of the membrane with an effective temperature 

f^ ' higher than that of the equilibrium membrane and an effective surface tension lower than that of 

oo ■ 

^^ . the equilibrium membrane. Asymmetric flip-flop leads to transient conformational changes of the 



membrane in the form of bud or blister formation, depending on the flip rate. 
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I. INTRODUCTION 

Biomembranes are non-equilibrium structures due to the non-thermal energy contribu- 
tions resulting from the activity of a wide variety of vicinal proteins. While the phase 
behavior and morphology of lipid bilayer membranes have been the subject of extensive 
amount of studies [l| , most of these studies have been on the equilibrium properties of the 
membranes. This has changed during the last decade or so with investigations of the effects, 
that active protein pumps have, on the undulations of lipid membranes, their morphology 
and the normalization of their mechanical constants [j, y, U, la, Ig, [Zl, |8|, |9| . Of particular 
importance is another class of membrane-bound proteins, that actively translocates phos- 
pholipids from one leaflet of a biomembrane to the other |lOl . I 111 . |l2j. The effects of these 
phospholipid translocators, on the mechanics and morphology of lipid membranes, has not 
received much attention. In this article, we present a study of self-assembled lipid bilayers in 
the presence of active lipid translocation, using dissipative particle dyanamics simulations. 

Lipid synthesis in eukaryotic cells takes place almost exclusively on the cytosolic leaflet 
of the endoplasmic reticulum (ER), which leads to an asymmetry in the lipid composition 
across the bilayer. In order to maintain a symmetric lipid density across the ER bilayer, 



nearly half of the newly synthesized lipids are rapidly translocated to the other leaflet |12 |. 
In contrast, the plasma membrane is marked by an acute asymmetry in the lipid compo- 
sition. Indeed while phosphatidylcholine and sphingomyelin are predominantly present in 
the exoplasmic leaflet, pho sphatidylserine and phosphatidylethanolamine are mainly found 
in the cytosolic leaflet [is!]. Maintenance of the symmetric lipid distribution in ER or the 
asymmetric lipid distribution in the plasma membrane cannot be mediated by thermally 
induced lipid movements (also termed passive flip-flops) alone. Indeed, passive flip-flops of 
phospholipids are energetically unfavorable due to the large energy barrier, 20-50 kcal/mol, 
associated with the translocation of the polar head group through the low dielectric per- 
mitivity hydrocarbon core of the bilayer. Consequently, the rate of passive flip-flops of 
phospholipids is extremely small, of the order of 10~^s~^ |l4l . Il5l |. i.e. on average, a single 
lipid experiences a thermally induced flip-flop every 24 hours. The lipid distribution across 
the bilayer is therefore actively jnaintained by a class of membrane-bound proteins known 



as phospholipid translocators [13]. These include the adenosine trh 
flippases and floppases, and the energy-independent scramblases [lO, ll 



losphate-dependent 



li 



Given the difficulties in purifying membrane-bound proteins, most of lipid transcolators 
have not been identified. Moreover, the mechanism(s) of active fiip-flop, mediated by phos- 
pholipid translocators, remain elusive. Few models have been proposed as mechanisms for 
active phospholipid translocations ly, [iTI] . In particular, Pomorski and Menon 12| recently 
proposed a mechanism similar to that of swiping a magnetic card through a card reader. In 
this model, the hydrophobic head group of the fiipped/fiopped lipid (magnetic strip of the 
card) is shielded from the hydrophobic environment of the bilayer, thereby facilitating its 
translocation across the bilayer. Recently, Sens ISj investigated theoretically the conforma- 
tional response of an infinitely large membrane to a localized disturbance in the form of a 
localized transbilayer asymmetry in the lipid density. There, he found that this asymmetry 
may transiently lead to the formation of a bud-like invagination followed by its relaxation. 
In this paper, we present a model for lipid translocation that is reminiscent of the magnetric 
swipe card model [l2|. The rnodel is then investigated via large scale dissipative particle 
dynamics simulations 19|, l20|, 



23|. To our knowledge, the presented work is the 
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first simulation study of the effect of active flip-flop on the mechanical and morphological 
properties of lipid membranes. 

II. BILAYER MODEL 

In the dissipative particle dynamics (DPD) model, for a self-assembled lipid bilayer in 
an explicit solvent used here, a lipid molecule is modeled as a flexible amphiphilic chain of 
beads consisting of one "head" (H) bead attached to three "tail" (T) beads via Hookean 
spring bonds. The solvent is modeled as single beads (W). All particles have the same mass 
m. In this model, interactions between any two non-bonded particles, within a range tq, are 
soft and repulsive. 

The forces acting on particles are grouped into three categories: (i) conservative forces, 
(ii) dissipative forces, and (iii) random forces. The conservative force between any two 
particles is 



p(C^i) 



C^ij ^ yij ) '■ij 1 



(1) 



where aij is the interaction strength between particles i and j, at respective positions r^ and 



j ' U ^ 3 ' ^-'^O. Yij 
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Vij\. Bonded particles belonging to a lipid also experience a 



conservative Hookean force given by 

F;j\ = -A;(l-r,,W&)rM+i, (2) 

where k is the spring constant and b is the preferred bond length. The dissipative force 
between particles i and j is given by 

FJf ^ = -'^ij^^ {rij ) {vij ■ Vij ) Tij , (3) 

where Tij is the dissipative strength for the pair {i,j), and Vjj = Vj — Vj is their relative 
velocity. The random force between i and j is given by 

Fjf = a.,(At)-V2^(r,,.)0,r.„ (4) 

where aij is the amplitude of the random noise for the pair {i,j), and Qj is a random 
variable with zero mean and unit variance which is uncorrelated for different pairs of particles 
and different time steps. Together, the dissipative and random forces act as a thermostat 
provided the fluctuation-dissipation theorem is satisfied. This yields to the following relation 
between T^j and aij, 

4 = 2TijkBT, (5) 

where k-Q is Boltzmann's constant and T is the thermostat temperature. In Eqs. ([T]), (^ 
and (j4]), the weight factor, uj{r), is chosen as 

1 — r/rn for r < rn 
coir) ^ ^ / - ° (6) 

for r > ro, 

where tq is the interactions cutoff. The particles trajectories are obtained by solving Hamil- 



ton's equations using the velocity- Verlet integrator 2J]. In the simulation, ro and m set 
the scales for length and mass, respectively. k-^T sets the energy scale. The time scale is 
given by r = [mrQ/k-QT] . The numerical value of the amplitude of the random force is 
considered to be the same for all pairs and is given by Oij = a = 3.0 (k'^T^m/rQ) , and 
the fluid density p = 3.0rJ"^. The amplitudes of the conservative force are chosen to be 
C(HH = ctTT = oiww = ctWH = '^'^k^T / tq and awr = olht = 2{){)k^T/rQ. In Eq. ([2]), the 
spring constant k = IOO^b^ and b = 0.45ro. The time step is chosen to be At = O.Olr. 
The fiat bilayer is initially constructed parallel to the xy - plane and placed in the mid- 
dle of the simulation box. It is then allowed to equilibrate until its normal fluctuations 



attain saturation. The total number of lipids used is 16,000 in a simulation box with size 
L X L X Lz = (86 X 86 X 40)rQ for the case of symmetric flip-flops and (80 x 80 x 46)ro for 
the case of asymmetric flip-flops. The system was subject to periodic boundary conditions 
in all three directions. 

III. FLIP-FLOP SCHEME 

We use the following two-step scheme for "flippase" action: (i) formation of a complex 
and (ii) translocation of a lipid from one leaflet to another. A lipid to be flipped is randomly 
selected from one of the two leaflets (refer to Figure [H for a schematic representation of 
the lipid to be flipped and the surrounding lipid molecules). Around this selcted lipid a 
fictitious cylinder is drawn which spans both leafiets of the bilayer. A fiippase complex is 
then defined as the set of all lipids inside this fictitious cylinder. The next step involves the 
action of a time- dependent fiipping force, F"(t) = F^(t)z, on the head bead of the selected 
lipid, so as to translocate it to the opposite leafiet. This force, F^{t), acts in the direction 
normal to the plane of the bilayer, and its magnitude is given by F^{t) = GAz{t), where 
Az{t) is the distance between the head bead of the lipid being translocated and the average 
2;-position of all the lipids, in the complex, in the opposite leafiet. This means that during 
the translocation process, the amplitude of the fiippase force decreases continuously with 
time. In order to conserve momentum within the fictitious cylinder, a force —F"-{t)/{Nc — l), 
where N^ is the number of lipid molecules in the fictitious complex, is concurrently applied 
onto the head beads of all other lipids within the fiippase complex. During the translocation 
process of the selected lipid, the head-tail repulsion of the selected lipid, with the other 
surrounding lipids in the membrane, is "screened" by temporarily choosing its amplitude 
anr = o^tt- This algorithm is therefore in-line with the recent swipe card model 12|. 



In Figure [2] the magnitude of the applied force, F^, normalized by G, is shown as a 
function of time, for the selected lipid as it translocates through the bilayer, mimicking the 
action of a fiippase. This figure depicts that the translocation time scale decreases as the 
amplitude, G, of the driving force is increased. In the remaining of the article, all results 
discussed are based on the case of G = lOk-QT/rQ, for which the typical time taken for 
fiipping a lipid from one layer is about r. During each time step, a number of lipids are 
selected at random. Flips are attempted with a probability, Pfiip, if the selected lipid is not 



already part of an active flipping complex. The success flip rate is measured by counting 
the number of lipids that reach the opposite leaflet in every time step. The flip and flop 
probabilities are respectively given by, 

^'^' ^ l + exp[-A{N^-N,-C)y ^^^ 



(8) 



and 

^^"■^ " l+exp[-A{N,-Nu + C)y 
where iV„ and A'^ are the numbers of lipid particles in the upper and lower leaflets, respec- 
tively. In Eq. ([8]), C controls the steady state number difference, AA^ = A^„ — A'^^^, and A 
is constant which fixes the width of the distribution of AA^. Having developed an active 
translocation algorithm, we now focus on the effect that active flip-flop has on flat mem- 
branes. We consider both cases of symmetric flip-flop, C = 0, and asynunetric flip-flop, 

IV. SYMMETRIC FLIP-FLOP 

Active flip-flop is labeled symmetric, if on average, the number of up-down translocations 
(flips) is equal to the number of down-up translocations (flops). At the beginning of the 
simulations, we start with a flat bilayer with exactly the same number of lipids in the two 
leaflets. We note that the rate of thermally induced flip-flops in this model is practically 
zero, in accord with experiments. Using the Irving-Kirkwood formalism [25], we calculated 
the lateral and normal components of the pressure tensor along the 2;-axis, and averaged over 
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271]. The tension of the membrane and its bending modulus are extracted 



the xy-plane 

from the structure factor of the out-of-plane fluctuations of the membrane height 201] 



A. Membrane Tension and Bending Rigidity 

The average orientation of the layer normal is taken to be along the z-axis. The steady- 
state profiles of the normal pressure Pn{z) and the lateral pressure Pl{z) along the bilayer 



normal are calculated from the pressure tensor using the Irving-Kirkwood formalism [25]. 
In Figure [H the normal and lateral pressure vs. z of equilibrium membranes are compared 
with that of membranes active symmetric flip-flop. The rate of flipping is 10 flips/ At over 



a membrane with projected area 86ro x 86ro. Although the details of the lateral pressure 
profile is dependent on the model used for the lipids, one observes that the fiip-fiop activity 
increases the lateral pressure in the two leaflets while the normal pressure is only weakly 
affected. As a result, active fiip-fiop reduces the effective tension on the membrane. 

By defining a height field h{x,y), which represents the position along the z-axis of the 
bilayer mid-plane at a point {x,y), and its Fourier transform /i(q) where q = {qxiQy), 
we calculate the circularly averaged structure factor S{q) = (|/i(q)p)/L^, where q = 
{qx + ly) ■ '^^^ lo^S wavelength deformations of a lipid membrane from its mean pla- 
nar conformation are well described by the Helfrich Hamiltonian [28|, 7Y[/i(x, y)] = 
f dxdy [f(V/i)^ + |(V^/i)^], where 7 is the membrane surface tension and k is its bend- 
ing modulus. The equipartition theorem of this model yields a structure factor. 

Hence, by plotting k-^T /q^S{q) as a function of g^, one extracts the tension on the membrane, 
7, from the intercept with the vertical axis and the bending modulus, k, from the slope at 
small wavevectors. This is shown in Figure H] for the cases of equilibrium and steady state 
fiip-fiop with varying fiip-fiop rates. Figure H] shows that as the rate of fiip-fiop is increased, 
the intercept of k-^T /q'^S{c[) is shifted to lower values, implying a reduction in the tension of 
the membrane. This is in line with the results from the lateral and normal pressure shown 
in Figure [31 However, Figure H] shows that the slope of k-^T /q^S{q) is independent of the 
fiip-fiop rate, implying that for the flip rates considered in our simulation, there does not 
seem to be any significant effect on the membrane bending modulus due to active fiip-fiop. 
We thus conclude that symmetric active fiip-fiops leads to an increase in the fluctuations 
of the membrane, which manifests itself in a decrease of the surface tension of the membrane. 

V. ASYMMETRIC FLIP-FLOP 

We now consider the effect of asymmetric fiip-fiop. In particular, we focus on the case 
where only lipids from the bottom leaflet are actively flipped to the upper leaflet. This 
is implemented by assigning a non-zero value for the constant C in the expression for the 
probability Pgip in Eq. ([8]). The value of C is kept large, equal to 10"^, such that the rate of 
accepted flips can be taken as a constant during simulation time. Flipping is restricted to a 



small square area (/ x /), termed the active region, in the central region of the membrane. 
This mimics the effect of flips due to localized flippases in a small region of a biomembrane. 
The effect of asymmetric flipping is then investigated as a function of / and the total number 
of flips per unit time (z/) . / is varied between lOro to SOrg and the flip rate is varied from 
u = 0.9 r-i to z/ = 20r-^ 

Asymmetric flip induces a finite difference in the lipid number densities, As = Su — Sd, 
where Su and Sd are the lipid number densities per unit of area, of the upper and lower 
leaflets, respectively. Furthermore, As increases with time as active flip proceeds. Therefore, 
asymmetric flip is characterized by an absence of steady state, and if maintained, leads to an 
instability of the membrane, in contrast to the case of symmetric active flip-flop. We found 
that asymmetric active flipping lead to the formation of two major transient morphologies 
corresponding to either buds or blisters, depending on size of the flip area and the flip rate. 

When the flip rate is low, u < 2r^^, we found that buds form. This is shown in Figure 
[5l In this case, a full bud is formed after about 800r. Once formed, the bud remains stable 
even after the flipping is stopped, and will relax back only in the time scale set by passive 
flips, which is much larger than the simulation time. The stability of the buds here has to 



do with the finite size of the membrane 
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30|. For an infinite membrane, the buds have a 
finite lifetime set by the relaxation of the in-plane density [l8|. 

When lipids are flipped at a relatively high rate, we observed the formation of blister 
structures in the active region. In the case with flipping rate z/ = 2.5r^^, it takes about 
300r for a blister to form. In Figure E], we depict snapshots showing the main stages of 
blister formation at this flip rate. The relatively high active flip rate and slow diffusion of 
the lipids results in high excess of lipids in the upper leaflet. This leads to the detachment 
of the upper leaflet from the lower leaflet, resulting in a protrusion that is reminiscent of 
a cylindrical micelle connected to the membrane. This is depicted in Figure [6]^a). The 
cylindrical micelle then grows into a sheet-like structure still connected to the membrane, as 
shown in Figure [6](b). If the process of active flipping is further continued, we found that 
the membrane becomes unstable. To avoid a destabilization of the membrane, we consider 
only the cases where active flipping is stopped once the sheet-like structure (Figure Et^b)) 
is formed. The relaxation of the blister depends on the stage at which active flipping is 
halted. If active flipping is stopped during the early stage of blister formation , shown in 
Figure [6](a), the blister relaxes back into the membrane. At late stages, when the sheet is 
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well formed (see Figure Mjo)) , in order to reduce the high edge energy, the blister curves 
forming a hook-like structure shown in Figure [6](c). Eventually, the sheet closes on to the 
membrane, forming a hemifusion state, reminiscent of that observed during the intermediate 



stages of vesicles fusion 3l|]. This is elucidated in the series of snapshots shown in Figure 
El^b-e). Eventually, the diaphragm in the hemifusion state ruptures into a pore as shown 
in Figure MS)- If active flip is stopped at an intermediate stage between those shown in 
Figure int^a) and (b), the blister relaxation is arrested at the hemifusion state shown in Figure 

VI. CONCLUSION 

The thickness of the bilayer obtained in the simulations is Avq. Comparing this with the 
thickness of lipid bilayers, 5 nm, we obtain rg ~ 1.25 nm. In phospholipid bilayers in the 
fluid phase, the diffusion coefficient of a lipid is typically D ~ lO^^^m^/s [l3|. Comparing 
this with the diffusion coefficient for lipids obtained in the current simulations, one estimates 
the DPD time unit, r ^ 0.2yus. For the case of G = S/c^T/rQ, we then find that the average 
duration of a single flip is about 2fis. Considering that the time scale for typical protein 
conformation change ranges from milliseconds to seconds, to make contact with real systems, 
one has to use smaller values of G. However, for practical reasons we use G = lOfc^T/rg. 
The typical simulation run in the present work is about 200/is corresponding to a diffusion 
length of approximately 15nm. Although the lipid density equilibrates very fast, the budding 
event is a result of the global asymmetry in the lipid number. 

In the case of an infinite membrane the formed buds will disappear on time scales larger 
than diffusion time. However in finite systems the resulting area difference of the two leaflets 
leads to buds that are stable within the time scale set by passive flip-flop rates |29|, |30| . 

The initial tension on the membrane also plays a role in determining the critical flip 
rate for bud and blister formation. In the case of membranes with flxed projected area, 
asymmetric flipping from the bottom leaflet to the upper leaflet leads to an asymmetry in 
the area per lipid in the two leaflets of the membrane. This in turn leads to an increase in 
the lateral tension of the bottom leaflet while that of the upper leaflet decreases. In order 
to equalize the area per lipid in the two leaflets, the membrane buckles, which will further 
increase the tension in both the leaves. Beyond a critical tension, the bottom leaflet will 



either rupture or decouple from the upper leaflet resulting in a blister. This also means that 
the flip rate at which blisters start to form should decrease with increasing tension of the 
initial equilibrated membrane. We confirm this in our simulations. 

In conclusion, we have presented a DPD model to study active flipping and its effects 
on a fluid bilayer membrane. We find that symmetric flip-flop results in a reduction in the 
tension of the membrane without much effect on its bending modulus. Asymmetric flip-flop 
results in non-equilibrium structures depending on the lipid flip rate. Slow flip rates involves 
membrane curvature and bud formation, whereas fast flip rates induces blister formation. 
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clearly visible in in (b), moves towards the cut plane in (b) to (e). The flip 
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FIG. 1: Flippase complex corresponding to the fictitious cylinder containing the lipid to be flipped. 
The lipid to be flipped is shown explicitly. However, only head beads of the other lipids in the 
flippase complex are shown. 
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FIG. 2: Normalized amplitude of the translocation force, F^{t)/G, vs. time for four values of the 
force magnitude, G. The data has been averaged over all translocated lipids and over time during 
steady state. 
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FIG. 3: Normal pressure Pn and lateral pressure Pl as a function of z (both equilibrium and with 
an attempted flip rate of 10 flips per At) for a system with L = 86ro. 
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FIG. 4: Structure factor vs. g^ for different values of attempted flip rates. Different symbols 
correspond to the cases of equilibrium (D), 1 flip/ At (■), 2 flips/ At (o), 5 flips/ At (•), and 10 
flips/At (A). The simulations are performed on a system containing 16000 lipids with L = 86ro 
and k-^T = 1.0. 
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FIG. 5: Bud Formation (only a small section of the membrane, with a cut through the budding 
region, is shown), a) Initial bending of both the leaflets and b) neck formation which eventually 
leads to a bud. The flip rate and size of the active region are v = 0.9r^^ and I = 10 respectively. 
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FIG. 6: Blister Formation (only a small section of the membrane, with a cut through the blister, 
is shown), a) Initial protrusion b) sheet formation c-e) fusion of sheet with membrane and f) pore 
formation. Note that the edge of the blister, clearly visible in in (b), moves towards the cut plane 
in (b) to (e). The flip rate and size of the active region are v = 2.5t~^ and / = 30 respectively. 
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FIG. 7: In its early stage blisters can also relax by folding into itself, forming a bulb-like structure. 
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